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A basic ingredient to understand large scale fluctuations in our present day Universe 
is the initial conditions. Using N-body simulations, we compare the clustering that 
arises from Gaussian and non-Gaussian initial conditions. The latter is motivated by a 
global texture model which has initial J-order correlations £j close to the strongly non- 
Gaussian regime ^ j ~ ^ 2 . The final amplitudes Sj = ^ j /t; 2 m the non-Gaussian 
(texture) model evolves slowly towards the (Gaussian) gravitational predictions but, 
even at ag = 1, are still significantly larger, showing a characteristic minimum with a 
sharp increase in Sj with increasing scales. This minimum, which is between 10 and 
15 Mpc, depending on the normalization, separates the regime where gravity starts 
dominating the evolution from the one in which the initial conditions are the dominant 
effect. In comparing this results with galaxy clustering observations, one has to take 
into account biasing, i.e. how galaxy fluctuations trace matter fluctuations. Although 
biasing could change the amplitudes, we show that the possible distortions to the shape 
of Sj are typically small. In contrast to the non-Gaussian (texture) predictions, we find 
no significant minimum or rise in Sj obtained from the APM Galaxy Survey. 



INTRODUCTION 

The large scale galaxy distribution can be used to study 
the origin and dynamics of cosmological fluctuations. The 
Q-Clust ering of matter density fluctuations 8r, smoothed over 
^ scale R, is characterized here in terms of the reduced J- 
5^ order moments fj(-R) = ^ <*> ^ . Assuming that gravity is 
the dominant dynamical effect, the evolution of £,j(R, t) is 
^completely fixed by the initial conditions ^j(K). For Gaus- 
sian inital conditions, i.e. {j = for J > 2, the leading order 
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gravitational evolution for small f 2 , 
-j- 



(R,t) 
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where f 2 follows the linear growth £ 2 (i?, t) = a(t) 2 ^ 2 (R) 
(Peebles 1980). Time independent gravitational amplitudes, 
Sj , have been predicted using perturbation theory (PT) and 
tested against N-body simulations (Fry 1984; Juszkiewicz et 
al. 1993; Bernardeau 1994; Baugh et al. 1995). Gaztanaga 
& Baugh 1995 have shown how the values of Sj change as 
a function of scale and time for different power spectrum, 
reproducing well the PT predictions on scales where f 2 < 1. 

Consider next the time evolution in the more general 
case of gravitational evolution from non-Gaussian initial 
conditions { j ^ 0. The leading contribution gives: 
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where Tj are geometrical factors. In an homogeneous dis- 
tribution — > in the limit f 2 — > 0, so that we can 



write -> (f 2 ) a , with a = a[J]. When a > J — 1 the 
initial conditions are forgotten, as the leading order effect 
of the evolution is the hierarchical term in equation (1). 
When J/2 < a < J — 1 we have quasi-Gaussian but non- 
hierarchical initial conditions. Evolution in { j have a domi- 
nant non-hierarchical term that grows as a , while the hier- 
archical term grows as a 2 ^ -1 - 1 and may not become signifi- 
cant until £ 2 ~ 1. Note that, as pointed out by Fry & Scher- 
rer 1994 that there is an additional non-Gaussian term that 
grows as a J+1 which for the skewness J = 3 contributes 
directly to S3 = { 2 /{ 2 so that S3 ^ S3 3 at all times. For 
J > 3 one expects a different scaling which could be used as 
a discriminating test (see Chodorowski & Bouchet 1996). If 
a < J/2 there are strongly non-Gaussian initial conditions 
that dominate the evolution as far as f 2 is small. 

The question we want to address in this letter (see Silk 
& Juszkiewicz 1991) is to what extend gravity is able to 
erase the trace of the initial conditions in the transition 
case to the strong non-Gaussian regime a ~ J/2. This is 
important because large scale galaxy surveys can be used 
to find such a trace (Gaztanaga 1992; Bouchet et al. 1993; 
Gaztanaga 1995) Previous work by Weinberg & Cole 1992 
provided interesting insights on different aspects of generic 
non-Gaussian models but did not address this question 
(other related works include Moscardini et al. 1993). Here 
we use a global texture model (Turok 1989) as a prototypical 
non-Gaussian model, which has recently been promoted as 
an interesting alternative for adiabatic inflationary structure 
formation (Cen et al. 1991). 
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Figure 1. The variance £ 2 a * two different epochs, <T8 = 0.1 
and <T8 = 1.0, for the texture simulations (closed triangles). Open 
circles correspond to the Gaussian CDM T = 0.5 at <T8 = 1.0. 
The dashed (continuous) lines show the linear PT predictions for 
the T = 0.5 (r = 0.7 ) CDM model at <r 8 = 0.1 and <r 8 = 1-0. 
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2 NON-GAUSSIAN TEXTURE MODEL 

The non-Gaussian global texture model dynamics is gov- 
erned by the evolution equation 



— + 2- — -V (*= -a , 

dr 2 a dr d<j> 



(3) 



where <f> is 4-dimensional scalar field, r is conformal time 
and V(<f>) = j(<f> 2 — <f>o) is the potential with the vacuum- 
energy value 00- We have used a non-linear sigma model 
for the texture dynamics, evolving fields with the modified 
leap-frog integrator algorithm (Pen et al. 1994; Nagasawa et 
al. 1993) and calculating the density perturbations at each 
time step. The full texture dynamics is stopped at as = 0.1, 
when we map all density perturbations by 100 3 particles 
to produce the initial conditions for a gravitational N-body 
simulation. The density fluctuations are then evolved by a 
P 3 M-code (Efstathiou & Eastwood 1981; Efstathiou et al. 
1988) until we reach as = 1.0. We have tested our simula- 
tions against the finite box size effects and our models are 
in good agreement with higher resolution simulations. The 
power spectrum P(k) for the texture model turns out to 
be similar to the adiabatic inflationary CDM model char- 
acterized by T = fi/t (Bond & Efstathiou 1984). Here we 
use Q = 1 and h = 0.5. Our simulations show that texture 
models can be scaled as adiabatic models simply by fi/t for 
different matter densities within 10% accuracy. More details 
are given in Mahonen & Efstathiou 1995. 



2.1 Clustering in the simulations 

The volume averaged correlations { j are estimated from mo- 
ments of counts-in-cells as described in Baugh et al. 1995. 
We average the results over 3 simulations in each ensem- 
ble and use the dispersion between members to estimate the 
sampling errors. Different stages in the evolution of the sim- 



Figure 2. The skewness S3 = C3/C2 111 the Gaussian T = 0.5 
CDM model (open circles) compared with the Non-Gaussian 
model (filled symbols). The top panel shows the initial condi- 
tions, 173 = 0.1, while the bottom shows time (73 = 0.4 (triangles) 
and (73 =1.0 (squares). Lines show the Zeldovich approximation 
(short-dashed) and the perturbation theory predictions for the 
r = 0.5 (continuous) and T = 0.7 (long-dashed) models. 



ulations are labeled by as , the linear PT variance in spheres 
of radius 8 A" 1 Mpc, i.e. al = a 2 J°(8). 

Figure 1 shows the variance in the texture model (filled 
triangles) at two different epochs, as = 0.1 and as = 1.0. 
Note that the initial results match roughly the linear vari- 
ance in the CDM T = 0.5 model (dashed-line), although in 
detail they are closer to the T = 0.7 shape (continuous line). 
The evolved non-linear variance at as = 1.0 is close to the 
corresponding non-linear variance in the Gaussian T = 0.5 
simulations (open circles). Note that the deviations from the 
linear growth (represented by the upper lines) are similar in 
both cases, indicating that the initial non-Gaussianities have 
only a small effect in the non-linear growth. 

The top panel in Figure 2 shows the initial values of 
the normalized skewness: S3 = { 3 /{ 2 in both models. The 
Gaussian model (open circles) matches well with the Zel- 
dovich approximation, as expected (see Baugh et al. 1995). 
The non-Gaussian model shows a characteristic increase of 
S3 with scale. At large scales, R > 10 ft -1 Mpc, a fit of the 
form f 3 = A^ 2 yields A ~ 1 and a ~ 3/2 + 0.1 (dotted line in 
Figure 2). This is close to the strongly non-Gaussian transi- 
tion mentioned in the introduction. The lower panel in Fig- 
ure 2 shows the evolution of S3 for Gaussian models (open 
circles) at as = 1 in comparison with the non-Gaussian mod- 
els at as = 0.4 (triangles) and as = 1 (squares). The Gaus- 
sian models reproduce quite well the PT predictions at large 
scales. As the texture simulations evolve, the shape and am- 
plitude of S3 in the non-Gaussian model slowly approaches 
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Figure 3. The kurtosis S4 = C 4 /C 2 as ' ln Figure 2. 
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Figure 4. Amplitudes Sj for J = 3 — 5 in the APM galaxies (sym- 
bols) compared with the matter amplitudes in the Non-Gaussian 
texture models normalized to <rg = 1 (short-dashed line), to 
as = 0.4 (dotted line) and scaled with a linear bias b = 1.5 
(continuous line) and non-linear bias b = 2.5 (long-dashed line). 



the one in the Gaussian models at small scales. At larger 
scales there is a change in the slope of Sj, showing a char- 
acteristic minimum which separates the regime where grav- 
ity starts dominating the evolution from the one in which 
the initial conditions are still the dominant effect, i.e. equa- 
tion (2). The scale at which the minimum occurs Rmm is a 
function of time (i.e. ag) moving from Rmm — 10 ft -1 Mpc 
at ag = 0.4 to Rmm — 15 ft -1 Mpc at ag = 1. 

A similar trend follows for higher order amplitudes Sj 
(see Figure 3 and Figure 4). We find that the initial condi- 
tions follow: f j = f 2 ^ + ' with e ~ 0.1 (dotted line in Figure 
3), again close to the transition to strongly non-Gaussian 
initial conditions. This tendency agrees well with the pre- 
dictions for J = 3 — 4 by Turok & Spergel 1991. The time 
evolution, at intermediate scales R ~ 10 ft -1 Mpc, can be 
approximated by: 

U*) = (l + s) - 2 ( J -%(* = 0) (4) 
between z = (ag = 1.0) and z = 1.5 (ag = 0.4). 

3 COMPARISON WITH THE APM 

In Figure 4 we show the values of S3 , S4 and S5 estimated 
from the angular APM Galaxy Survey (Gaztanaga 1994), 
assuming no evolution in Sj (closed symbols) or the tex- 
ture evolution (opened symbols) given by equation (4). The 
model for the evolution leads only to small differences, since 
the mean redshift in the APM is only z ~ 0.15. These three- 
dimensional estimates result from using a simple scaling law 
to model the projection effects. Although there are some 
potential problems with this modeling Bernardeau 1995, we 
believe that these results are accurate (see Gaztanaga 1995; 



Baugh & Gaztanaga 1996). Depending on the way counts- 
in-cells are estimated the mean angular amplitudes could 
increase or decrease rapidly with scale at the largest scales, 
i.e. I > 10° (> 40 ft -1 Mpc). These diverging estimates are 
not reliable as sampling effects from the finite APM volume 
dominate the statistics at these larger scales (Gaztanaga 
1994; Baugh & Gaztanaga 1996). 

Large scale galaxy fluctuations, 8 g , might be biased 
tracers of the underlying matter fluctuations, 8. To account 
for this possible bias and uncertainties in the normaliza- 
tion of the texture model, we consider different outputs and 
scale them with different biasing prescriptions. In Figure 
4 we show Sj in the the texture model for two different 
outputs, together with the values at ag = 0.4 normalized 
with a linear biasing relation, 8 g = b 8 which produces 
Sj ig = Sj/b J ~ 2 . We have chosen b = 1.5 as the optimal 
value to match the APM amplitudes Sj i9 around 8 ft -1 Mpc. 
Note that for ag = 0.4 the linear bias requires b = 2.5 if we 
want to fit the APM variance at 8 ft -1 Mpc, but this value 
of b produces a poor matching for Sj l9 . We introduce more 
biasing parameters with a non-linear transformation: 

8 g = f[8]^b[8 + ^8 2 + ^8 3 + ^8 i +...] (5) 

which for small variances, £ 2 < 1, still gives a linear relation 
for f 2 but changes the final amplitudes to S 9i j given by 
equation (10) in Fry & Gaztanaga 1993, e.g. S3 i9 = (S3 + 
3c 2 )/&. For the texture amplitudes at ag = 0.4 we have to 
fix b = 2.5, c2 ~ 1, C3 ~ 6 and C4 ~ —180 to match £ 2 , S3, 
S4 and S5 at 8 ft -1 Mpc in the APM. The resulting shapes 
are shown as long-dashed lines in Figure 4. For scales up 
to 40 ft -1 Mpc in S3, or up to 20 ft -1 Mpc in S4 we find no 
significant minimum or rise within the errors in the APM, 
in contrast to the unbiased or biased texture predictions. 
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4 CONCLUSION 

In Gaussian models with different initiai power spectrum 
there is an exceifent agreement for Sj between PT and N- 
body simulations on scales where the variance is approx- 
imately linear Juszkiewicz et al. 1993; Bernardeau 1994; 
Gaztanaga & Baugh 1995. The values of Sj do not evolve 
much with time. In contrast, we have seen here that the 
non-Gaussian (texture) model show a strong evolution from 
large initial values of Sj towards the values found in the 
corresponding Gaussian model (with similar initial power 
spectrum), although even at as = 1 there are important 
differences. The non-Gaussian models have a characteris- 
tic minimum that separates the regime where gravity starts 
dominating the evolution from the one in which the initial 
conditions are still the dominant effect. The scale at which 
the minimum occurs is a function of time (and therefore 
of the normalization) moving from Rmm — 10 ft -1 Mpc at 

(78 = 0.4 to Rmm — 15/t _1 MpC at (78 = 1. 

Our conclusions for the texture model seem different 
from those in Nagasawa et al. 1993, but a clear comparison 
can not be made as their results apply only to smaller scales 
and do not include the late-time gravitational dynamic that 
is simulated here with a full N-body simulation. 

A direct comparison of the mass amplitudes Sj with 
amplitudes in the APM, is not very useful, as the texture 
model does not even reproduce the second order statistics. 
The variance f 2 has less power at large scales than the APM 
galaxy distribution, which has a shape T ~ 0.2 (see Mad- 
dox et al. 1990; Gaztanaga 1995). A biasing between matter 
and galaxy fluctuations might account for part of this dif- 
ference at small scales, but at large scales biasing is unlikely 
to introduce distortions on the shape of f 2 . A texture model 
with more power on large scales (e.g. Q < 1), but with sim- 
ilar initial conditions, would produce similar results for 5j, 
but with slightly larger gravitational amplitudes, Sf . 

Biasing could also affect the final, directly observable, 
(galaxy) amplitudes 5j i9 (e.g. Fry & Gaztanaga 1993) and 
could even introduce scale dependence distortions for large 
biasing parameters (see Gaztanaga & Frieman 1994). Nev- 
ertheless, we have shown that these distortions are typically 
small in comparison with the pronounced non-Gaussian fea- 
tures that we find. Thus, one would still expect Sj to show a 
minimum or a steep rise around R ~ 15 ft -1 Mpc. For scales 
up 40 ft -1 Mpc in 53 or 20 ft -1 Mpc in 54, we do not find 
this characteristic minimum or steep rise in the APM fluc- 
tuations, in contrast to the predictions in the non-Gaussian 
(texture) model. 
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